Cauchy horizon stabihty in self-similar collapse: scalar radiation. 
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The stability of the Cauchy horizon in spherically symmetric self-similar collapse is studied by 
determining the flux of scalar radiation impinging on the horizon. This flux is found to be flnite. 

PACS numbers: 04.20.Dw, 04.20.Ex 



O 
O 
(N 

> 

o 
O 



> ■ 
in : 
m 
o 

o • 

^ ' 

o 

OH. 
I 



X 
S3 



I. INTRODUCTION 

Perhaps the richest source of examples of space-times 
admitting naked singularities is the class of spherically 
symmetric self-similar space-times. There is an extensive 
literature on the topic; the recent review of self-similarity 
in general relativity by Carr and Coley provides a suit- 
able bibliography. Of particular note in this class are the 
perfect fluid solutions studied by Ori and Piran the 
massless scalar field solutions studied by Christodoulou 
and by Brady |^ and the SU(2) sigma model solu- 
tions studied by Bizon and Wasscrman We mention 
these because (i) the matter model has particular inter- 
est, for either physical or mathematical reasons and (ii) 
these self-similar solutions are of interest in studies of 
critical phenomenon Q. More generally, self-similar so- 
lutions admitting naked singularities are of interest be- 
cause of what they may tell us about cosmic censorship. 
Intriguingly, the evidence is not all in one direction. Re- 
cent work has indicated the stability of perfect fluids ad- 
mitting naked singularities in the class of perfect fluid 
space-times Q], while for the case of the massless scalar 
field, generic spherical perturbations of self-similar initial 
data which correspond to naked singularities will lead 
to censored singularities Also, within the class of 
self-similar spherically symmetric space-times, the sec- 
tors corresponding to censored and to naked singularities 
are both topologically stable 

With these results in mind, the aim here is to begin a 
comprehensive study of the stability of Cauchy horizons 
in self-similar collapse. In the case of charged rotating 
black holes, the instability of the Cauchy (or inner) hori- 
zon has been firmly established (see for a review). 
This instability is in one way easily understood; an ob- 
server crossing the inner horizon views the entire history 
of the external universe in a finite amount of proper time, 
and so time-dependent perturbations of the exterior suf- 
fer an infinite blue-shift on crossing the horizon. This 
instability mechanism which can be "read off" the confor- 
mal diagram, does not have a counterpart in self-similar 
collapse which leads to globally naked singularities (see 
Figures 1 and 2). At best, one can speculate that the cur- 
vature at the regular center which diverges in the limit 



as the scaling origin is approached makes itself felt by 
perturbations approaching the Cauchy horizon. This is 
by no means convincing, and so a rigorous analysis is re- 
quired. We begin this analysis here by studying the prop- 
agation of scalar radiation in a fixed background (spher- 
ically symmetric, self-similar) space-time which admits a 
Cauchy horizon. 

In the following section, we define the class of space- 
times of interest and obtain some useful relations for the 
metric functions thereof. We consider spherically sym- 
metric space-times admitting a homothetic Killing vector 
field whose energy-momentum tensor obeys the dominant 
energy condition. (A complete account of energy condi- 
tions in spherical symmetry is given in the appendix.) 
For generality, no further restrictions are imposed at this 
stage, although some differentiability conditions at the 
past null cone of the scaling origin and at the Cauchy 
horizon will be imposed. Using co-ordinates adapted to 
the homothety and to the past null cones of the central 
world-line, simple conditions can be given on the metric 
which determine the visibility or otherwise of the sin- 
gularity at the scaling origin O. This allows a simple 
way of identifying both the past null cone M oi O and 
the Cauchy horizon TL. In Section 3, we determine the 
behaviour of completely general time-like geodesies (i) 
crossing M and (ii) crossing the Cauchy horizon. These 
are used to calculate fluxes of the scalar field at the re- 
spective surfaces. The minimally coupled scalar wave 
equation is studied in the next section. A mode decom- 
position relying on the Mellin transform is used, and the 
asymptotic behaviour of the general solution at M is de- 
termined. This is used to impose the boundary condition 
that an arbitrary observer with unit time-like tangent w° 
measures a finite flux v°'Va^\j^- We also demand that 
the influx at J~ be finite. The modes not ruled out by 
these boundary conditions are then allowed to evolve up 
to the Cauchy horizon and the flux w^Va^l^ is calcu- 
lated. Our principal result is that this flux is finite for 
all the cases we consider. 



II. SELF-SIMILAR SPHERICALLY 
SYMMETRIC SPACE-TIMES ADMITTING A 
NAKED SINGULARITY. 
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We will consider the class of space-times which have 
the following properties. Space-time {M.,g) is spheri- 
cally symmetric and admits a homothetic Killing vector 



2 




FIG. 1: A portion of the conformal diagram of the maximaUy 
extended Reissner-Nordstrom space-time. The observer O-jz 
falls through the event horizon (double line) and into the black 
hole. On crossing the Cauchy horizon (dashed) into a new 
asymptotically flat region, Or receives in finite time all the 
radiation emitted by C_e during its infinite history. 



r = 0, II > 




FIG. 2: Conformal diagram for an example of a self-similar 
space-time admitting a globally naked singularity. We use the 
advanced Bondi co-ordinates v and r described in Section 2. 
The Cauchy horizon is shown dashed, the event horizon as 
a double line and the apparent horizon as a bold curve. N 
is the past null cone of the scaling origin. Other structures 
can arise; there may be no apparent or event horizon; the 
censored portion of the singularity may be null; the naked 
portion of the singularity may be time-like. There is evidence 
that the naked singularity is generically globally naked. See 
for details. In every case for which the singularity is naked, 
the conformal diagram fails to display an obvious mechanism 
by which the Cauchy horizon may be destroyed, in contrast 
to the case illustrated in Figure 1. 



field. These symmetries pick out a scaling origin O on 
the central world-line r = (which we will refer to as the 
axis), where r is the radius function of the space-time. 
We assume regularity of the axis to the past of O and of 
the past null cone N of O. We will use advanced Bondi 
co-ordinates (w, r) where v labels the past null cones of 
r = and is taken to increase into the future. Transla- 
tion freedom in v allows us to situate the scaling origin 
at (u = 0, r = 0) and identifies v — i) with TV. The 
homothetic Killing field is 

£ - v—+r— 
dv dr 



The line element may be written 



^dn^ 



(1) 



where dil^ is the line element of the unit 2-sphere. 
The homothetic symmetry implies that F{v, r) = 
F{x),tp{v^r) = ^j{x) where x — v/r. The only co- 
ordinate freedom remaining in is v —>■ V{v); this is 
removed by taking v to measure proper time along the 
regular center r = 0. 

We will not specify the energy-momentum tensor of 
(A4,g), but will demand that it satisfies the dominant 
energy condition. A complete description of energy con- 
ditions in spherical symmetry is given in Appendix A. Of 
these, the following will be used. (These are equations 
(|A7|), (|A8|) and (|Al2p respectively.) 



xip' < 0, 
e'^iF' + xF^e'''^p') < 0, 
l-2F + 2x{F' + Fij') > 0. 



(2) 
(3) 
(4) 



We impose the following regularity conditions at the 
axis. As previously mentioned, we take v to be proper 
time along the axis for t; < 0. Noting that x ^ — oo on 
this portion of the axis, (|l|) then gives 



hm 2Fe^'f' = 1. 



(5) 



The other regularity condition that we use is that all 
curvature invariants are finite on r = 0, u < 0. In the 
present case, the (invariant) Misner-Sharp mass is given 
by 

E^^{1-2F). 

Then E/r^ is a curvature invariant; this term has the 
same units as e.g. R and ^'2. Demanding that E/r^ be 
finite on the axis yields 



lim F=-. 



(6) 



Combining (^ and (^ gives these regularity conditions: 



F(-oo) = -, V(-oo)=0. 



(7) 
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We define the interior region Mint of space-time to 
be the interior of M , i.e. the interior of the causal past 
of O. The exterior region Adext is defined to by = 
M.int^MiJM.ext- (These definitions are in line with those 
of Q.) We assume that the metric is regular throughout 
Mint U A/" - this set does not include O - by which we 
mean -0 e C^(— cx),0]. As a Cauchy horizon can only 
form in Mext, we assume further that Fjip ^ C'^ (—00, a;,) 
for some a;, > 0. As we will see, if a Cauchy horizon 
develops, it must be of the form x = Xc for some Xc > 0- 
Our assumption is that the metric is regular at least up 
until the Cauchy horizon. 

Since we are studying collapse, our assumptions must 
include some statement of regularity - in the sense of the 
absence of trapped surfaces - of an initial configuration. 
The 2-sphere (w, r) is trapped if and only if 

X{v,r) := g^^VarVbr < 0. 

In the present case, this is equivalent to -F < 0, and im- 
plies that the condition for an apparent horizon is _F = 0. 
So in order to express the notion that the matter is ini- 
tially in some non-extreme state, we rule out trapped or 
marginally trapped surfaces in the interior region Mint- 
We will also demand that Af is not foliated by marginally 
trapped surfaces, and so we take 

F(a;)>0 for aU x e {-oo,0]. 

Next, we point out the inevitability of there being 
a curvature singularity at O. Any curvature invariant 
which has units is of the form c(a;)r^^. For exam- 
ple. 



E 



1 - 2F 
2r2 ■ 



This term diverges as we approach O along the null line 
X — unless F{0) = i. But subject to the assumption 
that F > for a; < 0, we see that the surfaces x ^ Xc < 
are time-like. So we may also approach O along x ^ Xc < 
0, and we then see that E/r^ diverges unless F = i on 
{—00, 0]. Applying the same reasoning to the invariant 



2x{F' + F^')), 



regularity at O would require ^' = on (—00, 0] (we have 
used the boundary condition (j^) here). Hence Mint is a 
portion of flat space-time. So avoiding the trivial case 
implies the existence of a curvature singularity at O. 

Let us now prove the assertion above regarding the 
Cauchy horizon. 

Proposition 1 Let Xc be the first positive root of 
G{x) :— Fe"^' = if such exists and Xc — +00 other- 
wise. Then there are no future pointing outgoing radial 
null curves emanating from O in the region x G (0,Xc). 

Proof The outgoing radial null curves of (|l|) satisfy 

dr 



dv 



F{x)e^{x) = G{x). 



(8) 



Let {vi,ri) be a point on a solution curve 7^ of (|^) in 
the region < x < Xc- Then Xi = Vi/ri < Xc, and 
so G{xi) < ■ If a^c is finite, we note that a; = a;c is 
a solution of (^),and so by uniqueness, 7^ cannot cross 
X — Xc away from O i.e. for v > Q. Thus 



dx 
dv 



= -(1 -xG)>Q 
r 



for V £ (0, w,;]. Note that this inequality is immediate 
when Xc = +00. So the inequality applies generally and 
says that as w | 0, a; decreases and is bounded below by 
0. Hence the limit 



XI 



lim a;(w)|. 



< Xc 



exists and is non-negative. Thus either r — > > as 
f I - in which case the singularity is avoided - or r — > 
in the limit. In this case, 



xi 



lim — 

vio r 



— lim 



1 



i^io r'(v) 
lim ■ 



viQG{x) GixiY 

where all limits are taken along 7^ and I'Hopital's rule 
is used in the second line. The conclusion that xi < Xc 
is a root of xG = 1 contradicts minimality of Xc and 
completes the proof. □ 

Corollary 1 // G{x) < - for all x > 0, then the singu- 
larity O is censored. □ 

Corollary 2 // G{x) = i for some values of x > 0, then 
the Cauchy horizon Ti. of the space-time, where 
Xc is the smallest positive root of xG — 1. □ 

These results show an advantage of describing self- 
similar collapse in the co-ordinates v and r: the visibility 
of the singularity at O (and indeed the presence of an 
apparent horizon F = 0) can be read-off from the met- 
ric. More accurately, the presence of a naked singularity 
can be determined by tracking the evolution of metric 
functions, and without having to integrate geodesic equa- 
tions. 

An apparent horizon may form either before or after 
the Cauchy horizon. This horizon must be space-like, 
and the region lying to its future is trapped: 

Proposition 2 // F{xa) = for some Xa > 0, then 
X = Xa is space-like and the region x > Xa is trapped 
or marginally trapped. 

Proof Restricting to 'Sa ■ {v — Xar} in (|^) gives 



ds^ 



2a;ae'^(^")(l - XaGixa))dr'^ + r'^dn^, 
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which has spatial signature at G = F = when Xa > 0- 
From (^), we see that _F' < at an apparent horizon. 
Hence F{x) < for x > Xa- □ 
We conclude this section with a lemma which will play 
a central role in determining the stability of Ti with re- 
spect to scalar radiation. 

Lemma 1 G" < prior to the formation of a Cauchy 
horizon. 

Proof We note first that the results of Propositions 1 
and 2 show that < < 1 for x e {0,Xc). Then (||) 
gives 



and using (||) we get 

F' + F^' < F' + xF'^e^ip' < 0, 
i.e. G'{x) < for a; e (0,2;^). 
Corollary 3 G'(xc) < 0. 



□ 

□ 

We note that if G'{xc) = 0, then Rabk''k''\^ = 0, 
where is tangent to the outgoing radial null direction. 
This implies that there is no ingoing radiative flux of 
energy-momentum crossing the Cauchy horizon. We rule 
out this situation as being physically unrealistic and so 
we will assume that G'{xc) < 0. 



III. TIME-LIKE GEODESICS CROSSING Af 
AND n. 

The stability of the Cauchy horizon will be studied 
from the point of view of the behaviour of the flux of 
scalar radiation measured by an observer crossing the 
horizon. This flux is !F = w'^Va<&, where $ is the scalar 
field and is the unit tangent to an arbitrary time- 
like geodesic. Thus we will need to determine the be- 
haviour of the tangent w° for such arbitrary geodesies at 
the Cauchy horizon. Since we will impose boundary con- 
ditions on $ in terms of the fluxes at TV, we will need 
to do the same at this surface. The full set of equations 
governing time-like geodesies may be written in the form 



V - -{x{G' + G^') - ^')v^ - 

F' 



(9) 

2Ge'^w2 26"^^ + ^ = -1, (10) 
^ - ^, (11) 



where the overdot represents differentiation with respect 
to proper time r, L is the conserved angular momentum 
and f2 is an azimuthal angular variable, ( pi] ) plays no 
further role below, but is given for completeness. It is 
convenient to rewrite (|^) and (|l^) as a first order system. 



Defining X = {r,v,u)'^ where u 
may be written as 



i), these equations 



X = H(X) = 



^(2Ge^,2_|^_l) 
u 

l{x{G' + Gil;')-ij')u^ + ^e~^' 



(12) 

A future-pointing time-like geodesic crossing J\f corre- 
sponds to a solution of (p2|) with initial values > 0, 
Wo = 0, uo > 0. The assumptions of the previous 
section indicate that H is G^ in a neighbourhood of 
(rg, uo) £ R^, and so standard theorems imply the ex- 
istence of a G^ solution for X which exists for (at least) 
finite duration. Note that this implies that both v and r 
(via (p^)) are G^ functions of proper time r in a neigh- 
bourhood of J\f. Thus we can apply Taylor's theorem and 
write 

v{t) = MoT + tJaT^ +0(t^), 

r(r) — ro + riT + 0{t'^), 

where the coefficients V2 , ri can be given in terms of the 
initial data and metric functions and we have set r = 
at M. From this we may write down the following result 
which will be required below. 

Proposition 3 For any future-pointing time-like 
geodesic crossing M , we have 



V ^ UqT, V - Mo, 
^ Uo • ^ 

as T ^ where t — on the geodesic at J\f. 



(13) 
(14) 

□ 



Obtaining equivalent results at the Cauchy horizon is 
more difficult, as this corresponds to a singular point of 
the geodesic equations. Two things must be established: 
the existence of time-like geodesies crossing the horizon 
and the limiting values of the components of the tangent 
vector at the horizon. The proof below requires an as- 
sumption on the level of differentiability at the horizon 
which it would be desirable to remove. 

Proposition 4 Suppose that G and are differentiable 
at X — Xc- Then all radial time-like geodesies whose ini- 
tial points are sufficiently close to the Cauchy horizon 
will cross the horizon in finite time. For any time-like 
geodesic crossing the horizon, the components of the tan- 
gent X and V have finite non-zero values at the horizon 
which, denoting them by Xc and Vc respectively, satisfy 
the relation 



(15) 



2 Vcic 

where the subscript refers to the value of a quantity at 
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Proof (i) First, we establish a first order non- 
autonomous system for the geodesies. If is the homo- 
thetic KilUng vector field and u° is tangent to a time-like 
geodesic, then 



geodesic has initial point sufficiently close to x = Xc, 
cannot behave in this way. 

(iv) Consider a radial time-like geodesic for which v 
oo and x ^ Xc &s t —f oo. We have from 



where r is proper time along the geodesic (see e.g. Ap- 
pendix C of Q). Integrating yields 



■^^Cu'' = 2-e'^(l - xG)v - —e'f'x = k - t, 



9' 



for some k which is constant along the geodesic. Com- 
bining with 



-e'''(l-xG)v^ - ^e^vvx 
X x^ 



(16) 



(which is (10) written in terms of v and x) we obtain the 
first order system 



V 
1 X 



(17) 



V = --e-H^~xG)-\-{T-k)±Y^/% (18) 
2 V 



where 



2 r 2 

Y^{T-kf + 2— (1 - xG)e^(l -h ^). 
X 



We choose the upper sign, which corresponds to future- 
pointing geodesies. 

(ii) For radial [L = 0) time-like geodesies, (^ becomes 

v= (x^G' + xi;'{xG -1)) — . 

V 

Since G'{xc) < 0, the coefficient on the right hand side 
is negative for values oi x < Xc sufficiently close to Xc- 
Hence a geodesic with initial value xq = x(to) sufficiently 
close to Xc in this sense satisfies i; < for r > tq, and so 
V cannot diverge to infinity in finite time. 

(iii) Next, we establish that if w — > oo as r — > oo along 
a geodesic which does not cross the Cauchy horizon, then 
X Xc as r ^ oo. From (p^, we can write 

V X 

V x{l — xG) 

Integrating both sides yields and taking u = uq at some 
< < Xc, we get 



V <vo exp 



dy 



y(i - VG) 



for Xq < X. Thus if v diverges to -f oo, then so too must 
the integral. This can only occur if the integrand di- 
verges, i.e. if a; Xc- Now we show that provided a 



^^xlG'ixc) = -k^ <0 



as T — !■ oo. Integrating and reusing this relation yields 
the asymptotic relations 



V - civ^'' , (19) 
i k^civ^^''''-^ (20) 

oo for some ci > 0. Using these and (|l6|) we 



as r — 
obtain 



X 



C2V 



(21) 



as T ^ oo for some C2 > 0. We must have i as 
T ^ oo, for otherwise x is positive and bounded away 
from zero for an infinite amount of time and so x reaches 
Xc in finite time. Our present assumption is that this 
does not happen, so we must have fc^ < 1. 
The geodesic equations yield 



2Av + vBx + 2A'i 



vB'x^ = 0, 



(22) 



where 



A(x)^ -e^{l-xG), :re'^. 



Using the assumption that these terms are differentiable 
at the Cauchy horizon, we have from (|l|)-(|2|), fill 



2Ai 
2A'xv 
vB'x^ 



'2Ak^cW 



2k' 



^lu — o{v 

2A'ciC2V-^ = 0{v-^), 



Comparing these with (E4), we see that we must have 



We have 



lim vx 



'b 



V ^'-2 

iim v—x . 



-2x' 



Using the energy condition (g), we see that this term 
is strictly negative at the Cauchy horizon, and so this 
implies that x is positive for sufficiently large values of 
T. However this contradicts the fact that x(t) > with 
i ^ as T ^ 0. Hence the geodesic cannot extend 
to arbitrarily large values of v without first crossing the 
Cauchy horizon. 

(v) To conclude, it has been established that (at least 
some) radial geodesies cross the horizon in finite time 
and so with a finite value Vc of v. For any such geodesic, 
including non-radials, we can read off from (O) the non- 
zero and finite limiting value of x and from (|16[) we obtain 

dil). □ 
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IV. THE SCALAR FIELD ON THE CAUCHY 
HORIZON. 

Now we are in a position to examine the stability of 
the Cauchy horizon by measuring the flux of the scalar 
field in different regions of the space-time. 
In order to measure the flux J- = u°Va$ we need first 
the solution of the scalar wave equation, 



0. 



We exploit the spherical symmetry of the space time and 
split the scalar field, 

^v,x,9,^)=T{v,x)A{e,cf>), 

where we use the advanced null coordinate v, the ho- 
mothetic coordinate x, and the standard angular coordi- 
nates 0, (j). Then the line element in these co-ordinates 
reads 



= 2e-^[- 



G\dv' 



2e'f'v 



dvdx H — ^dfi^ 



By using separation of variables we arrive at a p.d.e. in 

V, X 



2x^ 



-2vT,^ -pe^T = 



(23) 



where p = 1(1 + 1) is the separation constant, I — 0,1,2 . . . 
is the multipole mode number, and ' denotes differentia- 
tion w.r.t. X. The complementary p.d.e. in 9, (p reduces 
to a form of Legendre's equation and is solved by the 
spherical harmonic functions, P"'-{9,(j)). 
We can perform a Mellin transformation on the field, de- 
fined by 



M{T}{x,n) = H^{x) 



T{v,x)v''-'dv 



which amounts to replacing T{v, x) with v"Hn{x), where 
n is an as yet unconstrained complex parameter. Equa- 
tion d23) thus reduces to an o.d.e. in H(x), 



2x' ( - 

X 



G]h" + {2n - 2x^G')H' 



2n 

X 



H = Q 



(24) 



where we have suppressed the subscript n. Performing 
the inverse Mellin transform on the solution of this o.d.e. 
over a contour in the viable range of n will return the 
solution to (p3|). 

This o.d.e. has a number of singular points, namely x = 
and the roots of xG = 1, the lowest of which we have 
defined to be Xc- The canonical form of a second order 
linear o.d.e. in a neighborhood ol x ^ xq is 

{x - xofH" + {x~ xo)q{x)H' + p{x)H = 0, 



and when we write equation (24) in its canonical form in 
the neighborhood of a; = 0, we find 



q{x) = 



n ~ x^G' 



1 



vG 



p{x) = 



2n + pe^x 
2{l-xG)' 



Since q{x) and p(x) are both G^ in a neighborhood of 
X = we can use the method of Frobenius to solve (24) on 
N ||l3| (see e.g. Chapter 3 of [|^). The indicial exponents 
are 1, —n. 

As it stands we cannot make any assumptions about n, 
however later analysis shows if — Re(n) > 1 the fiux of 
the scalar field will be always infinite on A/", thus we only 
consider — Re(rt) < 1. 

It is possible for 1 and —n to differ by an integer and so 
the method of Frobenius yields the following expression 
for the solution to (p4) in a neighborhood of a; = 0, 



m+1 



m=0 



-C2 < fclnx ^ ama;'"+^ + ^ 6„ 



(25) 



m=0 



m=0 



In this expression, ci and C2 are arbitrary constants, 
o-o — bf) = 1 with fc = if 1 and —n do not differ by 
an integer, ao = l,&o = with /c = 1 if 1 and ~n are 
equal, and — = 1 with k ^ Q ii 1 + n = m for some 
positive integer m. 

After some rearranging and some cancellations, the ex- 
pression for the flux on M is 



V 



n+m 



^ a„,{m + n+l)-^ 

m+n+1 



m=0 



V 



am{m + 1)- 



(26) 



J^2{v, r) 



V ^ 6m+i(m -I- 1)- 



m=0 
oo 

?n— 



bmim - n)- 



.m — n+l 



V 



ni—O 



Din 



OC 

-vk ^ 1^1-1- (to -I- n + 1) in (^-^ 



,rn— n+1 



,,n+m+l 



(27) 



where the 1 subscript denotes the ci part, and likewise 
the 2 subscript. 

The components of the velocity have been shown to be 
finite on J\f in Proposition 3, and we see that for the flux 
to have a flnite measure on J\f, that is when v = 0, we 
require 

Re(n) > 0. 

Under this condition we let the scalar field evolve 
towards Ti., and examine its flux there. 
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Note : A scalar field coming from past null infinity will 
have a finite flux thereon if Re(n) < 1. While this 
physically desirable condition should be imposed, it does 
not play any role in later analysis. 

When we write ( j^) in its canonical form around 
X — Xc, we find the coefficients are 




Now we reach an important distinction, whether G{x) 
has a unique lowest root or multiple lowest roots. We 
distinguish the two cases so: 

Lemma 2 When xG = 1 has a unique lowest root, 

xlG'{xc) + 1 > 0. 
When xG — 1 has a multiple lowest root, 

xlG'ixc) + 1 = 0. 

□ 

The two cases will lead to very different analyses, thus 
we treat them separately. 

(i) The first case leads to q{x),p{x) being G^ on 
X = Xc, thus Xc is a regular singular point and hence we 
can use the method of Frobenius. The indicial exponents 
are 0,1 ~ qa where 



xlG'{xc) - n 
xlG'{xc) + l' 



Since n > 0, Lemma 1 and 2 tell us qo < 0, hence 1 — go > 
0, which gives us 



H{x) ^GiY. Amr+^- 

^ oo oo 

+C2 hlnC ^ A„r+'"* + E 



90 



(28) 



m=0 



m=0 



where C = x — Xc, and the coefficients have the same 
structure as (pq). From this we calculate the flux. 



Ti = ^ (m + 1 - qo)AraC 

oo 



m=0 



^2 = if"C2 



/c^ An[lnC(™ + l-'Zo) + l]C 



m-qa 



'Vnv''-^G2 



oo oo 

^i3,„r + A; In c^A„,r+'"' 



m—0 



.(30) 



Using the flniteness of v, x given in Proposition 4, we see 
that if (Jo < 0, that is if n > 0, this expression is finite 
on Ti, i.e. when x — Xc — C — ^■ 

Thus in the case of xG — 1 having a unique lowest root, 
a scalar field measuring a finite flux entering the region 
will measure a finite flux on the Cauchy horizon. 

(ii) If x^G'{xc) + 1 = 0, Xc is an irregular singular 
point of (|^. Note that this is a special case which one 
would expect to correspond to a set of measure zero in 
the class of space-times under consideration. We label 
•q = Xc — X and examine solutions to the o.d.e. in the 
asymptotic limit 77 i (see e.g. Chapter 3 of @). 
We assume the solution to ( p^ ) can be written in the 
form 

reducing ( p4| ) to an o.d.e. in S. Now we assume the 
common property near irregular singular points, 

S = o{S^), ,7iO 

where the overdot denotes differentiation with respect to 
77. (|23) becomes a quadratic in S, 



{{xc -fl)- {xc - V?G} -(n+ {xc - vfc) 



S 



Xc-r] 2 

If we consider xG = 1 to have a lowest root of multiplicity 
fc, then we can write its Taylor series around 77 = as 

l-(xc- v)G(v) = P(v) = V*" ^ nf^k+i-, 

k\ 



This means if the lowest root is of multiplicity fc, we need 
the metric functions to be G^. This is not too much of 
a restriction however, since the class of functions with 
roots of multiplicity k becomes very small as k increases, 
meaning we are dealing with a very special case in this 
analysis. 

We can make the approximation 

n+{xc-vfG ~ n + 1, i 0, 

and since we assume the metric coefficients are at least 
, we can approximate e"^ by the first term in its expan- 
sion, Co, in the limit 77 | 0. Thus we arrive at a quadratic 
in 5, 



Tj'iSy-aS^(3, 77 J, 0, 



(32) 



k\{n + l) 

XeP«(0)' 



k\ f n pcQ^ 

XcPW(O) [~c^~, 
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where a, /? > (if Re{n) > 0) and constant in the Hmit 
r] I 0, and k > 1. This quadratic has two solutions 
corresponding to two Unearly independent solutions of 
(^), which are 



Si 



(fc-1) 



'a"^ ' (fc + 1) 
At this point we verify our earlier assumption, namely 

Thus we have constructed two solutions to (E4), 



Hiiv) 



r/^ exp 



(fc-1) 



exp{0(?7)} 



(33) 
(34) 



Both of these functions and their derivatives are finite 
in the limit ?7 J, 0, a; — > a:c if Re{n) > 0, and thus the 
resulting expressions for the flux arc finite, where again 
we use Proposition 4. 
We summarize thus; 

Proposition 5 Let space-time (A4,g) satisfy the require- 
ments o/ Section II and admit a Cauchy horizon x — Xc- 
Assume also that gab € at x = Xc- Then a scalar field 
which has a finite flux on M , the past null cone of O, will 
also have a finite flux on the Cauchy horizon, TL. 



V. CONCLUSIONS 

We have shown that the Cauchy horizon formed by col- 
lapse in a self similar, spherically symmetric space-time 
is stable with respect to scalar radiation. This space-time 
is very general, the only other constraints being that the 
field satisfies the dominant energy condition, and, other 
than the special case discussed in Section IV(ii), we re- 
quire the metric functions to be on M and TL. These 
differentiability conditions are stronger than one would 
like to assume (c/. the C° Cauchy horizons appearing in 
the collapse of wave maps Q), but are as low as one can 
go without having to resort to a generalised solution con- 
cept for the wave equation. 

The next step is to examine whether linear perturba- 
tions of the metric functions will lead to an unstable 
Cauchy horizon, as is seen, for example, in the Reissner- 
Nordstrom solution. Such an examination would be more 
significant in considering cosmic censorship. Is it difficult 
to anticipate the general outcome of such an examina- 
tion. One expects to observe instability for the case of 
a massless scalar field [||, but stability for (some sectors 
of) perfect fluid collapse 0. The present results and the 
Cauchy horizon stability conjecture would lead one to 
expect stability in general p5| . 
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APPENDIX A: ENERGY CONDITIONS IN 
SPHERICAL SYMMETRY 

1. Spherical symmetry 

We write the line element in double null coordinates; 

ds^ ^ -2e-^f dudv + r^dn'^, 

where f — f{u,v), r — r{u, v) and dfl^ is the line element 
on the unit 2-sphere. The non-vanishing Ricci tensor 
terms are 



Ruu 


— 2r (t^uu 


h 2r,/„), 


R-vv 


= -2r-i(r™ 4 


- 2r,/,), 


Ruv 


— — 2r {^T"uv ~ 




Ree 


= csc^ 6R^^ = 


2 2c ^fTiiy^ 



where E is the Misner- Sharp mass, 

E ^'^{l + 2e^fr^r,). 

Subscripts on /, r denote partial derivatives. The only 
non- vanishing Weyl tensor term is 

*2 = -~ + \e^^\fuv+r-\uv) 

6 7"' 3 



E 1 



where x^ are co-ordinates in the Lorentzian 2-space. 



2. The strong energy condition 

Our aim is to write down a set of conditions on the 
curvature terms listed above which are equivalent to the 
strong energy condition: 

RabV^v'' > 

for all (future-pointing) causal vectors v. 



a. Null vectors 

The radial null directions are <5°, 5". These give 
Ruu > 0, Rvv > 0. 

At any point, the general non-radial null direction may 
be written as 



(Al) 
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The null condition is 



Time-like vectors 



a(3 = 2r'^e'^f-/^. 



We find 



E 



RabV^v'' = a^Ruu + /3^Rvv + 2a/3(2/™ + 2-36-2/). 

This is non-negative for all non-radial null vectors if and 
only if it is non-negative for all values of a, (3 with a(3 > 0. 
In turn, this is true if 

min (5(/i) > 0, 

Ai>0 



where fi = a/ f3 > and 



E 



Again we write 

and we can use the time-like condition gabV°'v^ = —1, so 
that 

72 ^ r-2(2e-2/a/3- 1) > 0. 

So in this case we are minimizing over the set afi > ^e^^. 
we do this by minimizing over the hyperbola a/? = c 
and then minimizing over all hyperbolas c > ■ This 
yields the conditions above and the extra condition 

'^{_RuuRvv) ^ + Juv f T"uv ^ 0. (-^5) 



If Ruu — Rvv = 0, this is simply 

If Ruu — and Rw 7^ 0, the condition is equivalent to 

E 

2^{2Uv + 2^e-2/) + R^^>o 

for all > 0. This is satisfied if /„„ + ^e'^/ > 0. If 
fuv + ^£-2/ < 0, then the condition will be violated 
for sufficiently large values of /i (which can always be 
chosen). The same holds for Ruu ^ 0,Rvv = 0. Thus if 
RuuRvv = 0, the strong energy condition holds for null 
directions if and only if 



fuv + ^e-^f > 0. 



So now assume that Ruu > 0, Rw > 0. The quadratic 
Qip) has a global minimum at 



f^* - -'^RuuUu 



E 



while Q(0) = R„, > 0. Thus Q(/x) > for > if and 
only if either ^, < or ^* > and Q{fJ,*) > 0. 

/Xh. < if and only if /„„ -I- p-e^^/ > Q. 

/X* > if and only if /„„ -I- -^-e-^/ < Q. In this case, 



^ \fuv + < ^(^u«i?™)'/' 

^ f + —e-^f > --(R R )i/2 
r-' 2 

Combining these results, we can say: 
RabV°'v^ > for all null w° if and only if 

Ruu > (A2) 
Ruu > (A3) 

l-iRuuRvv)^^^ + fuv + ^e-^f > 0. (A4) 



3. The weak energy condition 

The weak energy condition TabV°'v^ > for all 
causal w° can be written, using Einstein's equation, as 
RabV^'v^ > eR/2, where e — gabv°'v''. {R = Ricci scalar.) 
Thus the only extra work to do is for time-like vectors. 
The algebra involved in the previous section only needs 
minute changes, and we can show that the weak energy 
condition is equivalent to (A2) - (A4) and 



2 y^uu^vv } 



%e-^f > 0. 



(A6) 



4. The dominant energy condition 



This states that for every future-pointing timelike v"", 
the vector —T°'^vi, is non-spacelike and future-pointing, 
and Tabv'^v'' > 0. Using the usual general form for v"", 
we again have quadratic inequalities for the parameters 
a and /3 which can be treated in the usual way. (The 
non-spacelike condition is gacT°''^T"^VbVd < 0; the left 
hand side is homogeneous of degree 2 in (a, P) and so 
quadratic in fi.) The future-pointing condition is simple 
to examine by assuming that u, v increase into the future. 
The resulting inequalities are 









> 


0, 






Rvv 


> 


0, 


E 


2/ _ 


^ fuv 


> 


0, 


E 

13^ 


2/ _ 


f fuv 


> 


1 fuv ~ f fuv 



Using the first three of these, we see that the left hand 
side of the fourth is non-negative, and so the fourth is 
equivalent to the two inequalities 

-(R R )^/^ + —e~^^ + 2r~^r -f >0 

i-. V ^uu-'-f'VV J ^ „ c ^1 I uv J uv _ 

2 r-^ 

(-^uii-^uv ) ^ ~l~ ^ ~\~ fuv ^ 0- 
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Note how (as expected) some of these are the same as 
some of the strong and weak energy conditions. 



Summary: covariant form of the energy 
conditions 



The energy conditions are given here in terms that use 
r^v and fuv A more transparently covariant form is 
obtained by using R and ^"2 instead of these two. Then 
the resuhs are as follows (we note that the signs of 
Ryy and the term e^^ RuuRw are invariants - the last of 
these is defined in terms of contractions of Ricci with the 
two invariantly defined radial null directions): 

The strong energy condition is equivalent to 



Ru 
Rv 



> 0, 

> 0, 



(A7) 
(A8) 



2 r-^ 12 

ie2/|i?„„i?,„|i/2 + 4 + *2-^ > 0. (AID) 
2 r-^ 



The weak energy condition is equivalent to ( A7) 
(^) and 



L^f\R^^R^^\^/^ + ^ + ^,^ + ^ > 0. (All) 
2 r-^ 12 



The dominant energy condition is equivalent to 
(M), (pi) and 



^ + *2 + ^>0, (A12) 
ie2/|i?„„i?,„|i/2 + | > 0. (A13) 
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